First-principles study of the structural energetics of PdTi and PtTi 
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The structural energetics of PdTi and PtTi have been studied using first-principles density- 
functional theory with pseudopotentials and a plane-wave basis. We predict that in both materials, 
the experimentally reported orthorhombic _B19 phase will undergo a low-temperature phase transi- 
tion to a monoclinic _B19' ground state. Within a soft-mode framework, we relate the -B19 structure 
to the cubic B2 structure, observed at high temperature, and the 519' structure to _B19 via phonon 
modes strongly coupled to strain. In contrast to NiTi, the -B19 structure is extremely close to hep. 
We draw on the analogy to the bcc-hcp transition to suggest likely transition mechanisms in the 
present case. 
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I. INTRODUCTION 

Shape memory alloys (SMA) have attracted a great 
deal of attention due to their important technological 
applications, including mechanical actuator devices and 
medical stents. The shape memory effect also gives rise 
to superelasticity, which finds applications in stents and 
spectacle frames. 

The shape memory effect is related to a reversible 
martensitic (diffusionless) phase transformation. It has 
been shovi^n that the martensitic transformation can be 
induced by applied fields, temperature or both, and the 
mechanical properties of materials, therefore, can be con- 
trolled accordingly. In many systems, including those 
discussed in the present work, alloying can dramatically 
change the properties and transition temperatures of the 
materials, reflecting the importance of electronic fea- 
tures, specifically Fermi surface effects, in the structural 
energetics of SMA. 

There are several complementary approaches to mod- 
elling of the shape memory effect. Continuum mod- 
elling allows investigation of the microstructural behav- 
ior, specifically martensitic twins, at the relevant long 
length scales. Material-specific behavior is incorporated 
through an empirical functional for the free energy in 
terms of strain and a twin boundary energy to set the 
length scale. In atomistic models, the individual atoms 
are considered explicitly and their interaction given by 
an interatomic potential, which may be determined em- 
pirically, from first-principles density-functional-theory 
(DFT) calculations, or a combination of the two. Crys- 
tal symmetry and defect energies emerge from this ap- 
proach, which gives microstructures with both natural 
length scales (from defect energies) and time scales (since 
the atoms have definite Newtonian forces and masses). 
However, in atomistic models, the electronic degrees of 
freedom do not appear explicitly. First principles DFT 
methods are so computationally intensive that direct 
studies of microstructural behavior are impossible, but 
they are valuable both for obtaining quantitative atomic- 
level information regarding energies, forces and stresses 
independent of empirical input, and for understand- 



ing the electronic origin of this behavior. Thus, first- 
principles investigation of the energetic instability of the 
high-temperature structure towards the low-symmetry 
martensitic structure is in itself quite illuminating. The 
resulting information can then also be used as inputs to 
atomistic [Q and continuum modelling of shape memory 
behavior. 

Typically, martensitic transformations are described 
using the strain as an order parameter, the classic exam- 
ple being the Bain bcc-fcc transformation of iron. How- 
ever, there is an alternative approach appropriate for 
cases where the strain degrees of freedom are coupled to 
atomic degrees of freedom (phonons). Following the soft- 
mode theory of structural transitions, Q we start from a 
high-symmetry reference structure (here B2) and freeze 
in unstable phonons of this structure, with corresponding 
lattice relaxation, to produce the ground-state structure. 
The symmetry of the phonons determines the symme- 
try of the low temperature structure. This approach has 
been successfully used in the study of minerals^, ^, ||] 
and ferroelectric materials |||, ^, ^ and has been ex- 
tended to shape memory alloys in our previous study of 
NiTi0. 

Closely related to NiTi, PdTi and PtTi are shape 
memory materials with reasonable structural simplic- 
ity and extraordinary mechanical behavior. They un- 
dergo a martensitic transformation at tunable temper- 
atures: PdTi transforms at 810K, but this can be re- 
duced to 410K with 8% substitution of Cr for Pd0. 
The high-temperature "austenitic" phase has a simple 
cubic B2 structure (space group PmSfh), while the am- 
bient temperature "martensitic" phase has been reported 
as the orthorhombic B19 structure |l3| (space group: 
Pmma). Previous first-principles studies in PdTi and 
PtTi|14| 1^, have shown that the observed electronic 
and elastic properties of the B19 structure are well repro- 
duced by density- functional theory calculations assuming 
the experimentally determined structure. 

In this paper, we investigate the structural energet- 
ics of PdTi and PtTi from first-principles calculations of 
phonon frequencies as well as total energies. This allows 
us to examine local as well as global stability and to in- 
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vestigate transition mechanisms, drawing on the strong 
analogy between the B2-B19 and bcc-hcp transforma- 
tions and showing that couphng of unstable modes to 
the strain is a crucial aspect of the structural energetics. 
In Sec. II, we describe the first-principles calculations. 
In Sec. Ill, we present and discuss the results for the 
phonon dispersion of PdTi and PtTi in the B2 structure 
and for the relaxed structures in which unstable modes 
and strains are coupled, yielding a low-symmetry ground 
state. In addition, we present results of calculations of 
the electronic structure, identifying and discussing fea- 
tures that affect the relative stability of the phases. Sec. 
V concludes the paper. 



II. CALCULATIONS 

First-principles total energy calculations were carried 
out within density-functional theory with a plane-wave 
pseudopotential approach. The calculations were per- 
formed with the Vienna ab-initio Simulations Package [p^ 
p^ , using the Perdew-Zunger|[T9| parametrization of the 
local-density approximation (LDA). Vanderbilt ultrasoft 
pseudopotentials pof were used. Our pseudopotentials 
include nonlinear core corrections and for Ti, we treated 
the occupied 3p levels as valence. The electronic wave 
functions were represented in a plane-wave basis set with 
a kinetic energy cutoff of 278eV. The Brillouin zone (BZ) 
integrations were carried out by the Hermite-Gaussian 
smearing technique pl[ | with the smearing parameter of 
O.leV. The unit cells contain two atoms in the cubic B2 
structure and four atoms in the orthorhombic B19 and 
monoclinic B19' structures. The calculations were per- 
formed with al6x 16x16 Monkhorst-Pack (MP) fc— point 
mesh for the cubic B2 structure and a 12 x 12 x 16 MP 
fc— point mesh for both orthorhombic B19 and monoclinic 
B19' structures (space group: P2i/m), corresponding 
to 120 k points in the ^ irreducible BZ of the sim- 
ple cubic cell, 288 k points in the | irreducible BZ of 
the orthorhombic cell and 576 k points in the j irre- 
ducible BZ of the monoclinic cell. This choice of param- 
eters converges the total energy to within 1 meV/atom. 
The density of states (DOS) for the B19 and B19' struc- 
tures were calculated using the tetrahedron method with 
Blochl corrections!^. The BZ's for the orthorhombic 
B19 and monoclinic B\9' structures are different. To 
compare the band structure of the two structures, we la- 
bel the B19 band structure by regarding the B19 struc- 
ture as a special case of B19' . 

The phonon dispersion relations were obtained with 
the PWSCF and PHONON codes ||, using the Perdew- 
Zunger|Tg[| parametrization of the LDA, as above. Ul- 
trasoft pseudopotentials ||2^ for Pd, Pt and Ti were 
generated according to a modified Rappe-Rabe-Kaxiras- 
Joannopoulos (RRKJ) scheme ||2^ with three Bessel 
functions |25| . The electronic wave functions were repre- 
sented in a plane- wave basis set with a kinetic energy cut- 
off of 408eV. The augmentation charges were expanded 



up to 9000eV. The Brillouin zone (BZ) integrations were 
carried out by the Hermite-Gaussian smearing technique 
PH using a 56 fc-point mesh (corresponding to 12 x 12 x 12 
regular divisions along the kx, ky and k^ axes) in the ^ 
irreducible wedge. The value of the smearing parameter 
was (T=0.2eV. These parameters yield phonon frequen- 
cies converged within 5 cm^^. The dynamical matrix 
was computed on a 6 x 6 x 6 g-point mesh commensurate 
with the fc-point mesh. The complete phonon dispersion 
relation was obtained through the computation of real- 
space interatomic force constants within the correspond- 
ing boxlUl. 

The choice to use two different first-principles codes 
was dictated by the individual strengths of each. VASP 
has a highly efficient scheme for calculating total energies, 
forces, and stresses, and relaxing to the minimum en- 
ergy structure, but does not have the density-functional 
perturbation theory capabifities of PWSCF/PHONON. 
Even with slightly different pseudopotentials and A;-point 
sampling, the results of the two codes are quite compat- 
ible. For example, the difference between the computed 
lattice parameters for the B2 structure of PdTi is less 
than 0.2%, and for PtTi the difference is less than 0.1%. 
Comparisons of normalized eigenvector components com- 
puted by VASP using the frozen phonon method and by 
PWSCF/PHONON also show good agreement, generally 
within 5%. 



III. RESULTS 

For the cubic B2 structure, our calculations yield the 
equilibrium lattice parameters of 3.112A and 3.125A 
for PdTi and PtTi respectively. For comparison, 
we also performed full-potential linearized-augmented- 
plane-wave calculations (FLAPW) within the LDA[p7[. 
The results are in excellent agreement with FLAPW re- 
sults of 3.113A (PdTi) and 3.127A (PtTi) and in in good 
agreement with experimental values of 3.18A and 3.17A. 

The phonon dispersion relations along high symmetry 
lines, computed at the theoretical lattice parameters, are 
shown in Figure |l|. The frequencies are obtained by tak- 
ing the square root of the eigenvalues of the dynamical 
matrix p9|. Imaginary frequencies, as plotted as negative, 
are obtained from negative eigenvalues of the dynami- 
cal matrix. Thus, the structure is dynamically unstable 
against distortions following the corresponding eigenvec- 
tor. It should be noted that the phonon frequency lj is 
not the reciprocal of the period of oscillation of this mode 
(as measured in molecular dynamics) nor is it the energy 
difference between adjacent quantum levels (as measured 
in neutron scattering experiments). These three quanti- 
ties are equal only for a stable harmonic crystal. In the 
materials studied here the unstable modes may be related 
to a soft mode||] (as defined by MD or neutron scatter- 
ing) at high temperature, and even the stable modes are 
expected to be significantly renormalized as a function of 
temperature by anharmonic effects. 
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FIG. 1: Phonon dispersion relations for PdTi (upper) and 
PtTi (lower) in the B2 structure calculated at the LDA equi- 
librium parameters 3. 112A and 3. 125A respectively. The neg- 
ative slope of the acoustic V — M branch corresponds to a 
pure elastic instability (c' = ^(cii — C12) < 0). Symmetry 
labels are assigned according to the conventions of Ref. Eq 
with Pd/Pt at the origin. The imaginary frequencies of the 
unstable modes are plotted as negative values. 



The dynamical matrices are related by mass factors 
to the force constant matrix: the second derivatives of 
the internal energy with respect to atomic displacements. 
The eigenmodes of the force constant matrix describe the 
potential energy landscape, and a negative eigenvalue in- 
dicates a static instability against a distortion following 
the corresponding eigenvector. While the actual normal- 
ized displacements of these eigenmodes are in general 
slightly different they carry the same symmetry labels as 
the eigenmodes of the dynamical matrix. Either choice 
is expected to serve as a useful pointer to a lower energy 
structure if the distorted structure obtained by "freez- 
ing in" an unstable mode is relaxed using first-principles 
forces and stresses, as we describe below. 

The phonon dispersion relations shown in Figure |l| 
show instability of the B1 structure similar to and even 
stronger than that of NiTi[|l0|. There are large regions 
of reciprocal space where one, two or even three modes 
are unstable, with dominant instabilities at M and along 
T-i?. The phonon instability shows that the observed 
high-temperature Bl phases of PdTi and PtTi are dy- 
namically stabilized by anharmonic phonons, and should 



be characterized by large fluctuating local distortions. 
The calculated phonon dispersions are also reminiscent 
of those of unstable bcc materials such as Zr and Ti, 
which undergo martensitic transformations to hep or 
a;[3C|] (via the M and T-i? bcc-phonon equivalents re- 
spectively) phases. The analogy based on the view of 
519 phases of PdTi and PtTi as chemically ordered hep 
will be further strengthened below. 

In the soft mode approach, we search for local energy 
minima by choosing an unstable mode of the high sym- 
metry structure, freezing in the distortion with varying 
amplitude, and relaxing the resulting structure. In many 
cases, the mode with the largest negative eigenvalue will 
generate the lowest energy structure. However, this is by 
no means generally true, as the energy gain is determined 
not only by the curvature of the energy surface but by 
higher order terms as well as the strength of coupling 
to strain and other modes, both unstable and stable, of 
appropriate symmetry. Indeed, in PdTi this "most un- 
stable" mode (i.e. largest negative eigenvalue) lies in the 
T-i? branch rather than at M. This mode is typically un- 
stable in dynamically-stabilized bcc materials such as ti- 
tanium and zirconium where it is associated with a phase 
transition to the ideal (C32) phase. The fact that in 
the chemically ordered analog (Pd/Pt)Ti, the observed 
phase transition is to i319 rather than C32 illustrates the 
importance of anharmonic effects and strain coupling in 
the energetics of these materials. 

In the Bl structure, a doubly-degenerate unstable Mg 
zone-boundary mode implies a doubling of the unit cell 
to a=(001)hcc, b=(110)bcc, c=(lTO)bcc- For the distor- 
tion to be frozen in, we choose the eigenvector that gives 
the space group Pmma of the B19 structure (Figure ||). 
The unit cell is orthorhombic, with lattice parameters a, b 
and c, two Wyckoff positions occupied by Ti and Pd/Pt 
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FIG. 2: Eigenmode of the doubly-degenerate M5 unstable 
phonon in the B2 structure. This mode generates the 519 
structure in PdTi and PtTi. The structure is shown projected 
along the z direction, with Pd/Pt represented by filled circles 
and Ti by open circles. 
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TABLE I; Computed structural parameters and total energies of PdTi and PtTi in the B19 and 519' structures, compared 
with the computed Llo structure, the computed bco structure, the experimental _B19 structure, with previous calculations]!^] 
and with two special _B19 structures with b/a and c/a corresponding to ideal hep and B2 (bcc) structures. Volume is given in 
A^/formula unit, and energy in eV/atom. Wyckoff positions refer to space group P2i/m (_B19'). 





Structure 


Volume 


h/a 


c/a 


7(°) 


Wyckoff position 


X 


y 


z 


E — Ebi 


PdTi 


519 


29.63 


1.753 


1.634 


90 


Pd(2e) 
Ti(2e) 




0.5 


0.6866 
0.2008 


0.25 
0.25 


-0.0917 




B\Q' 


29.64 


1.758 


1.633 


93.39 


Pd(2e) 
Ti(2e) 


0.0114 
0.4475 


0.6827 
0.1964 


0.25 
0.25 


-0.0924 




Llo 


29.32 


1.000 


1.375 


90 










-0.0870 




bco 


30.16 


1.837 


1.579 


105.79 


Pd(2e) 
Ti(2e) 


0.0882 
0.3491 


0.6764 
0.1983 


0.25 
0.25 


-0.0573 




B19 " 


31.33 


1.74 


1.62 


90 














B19 


31.74 


1.75 


1.64 


90 














B19 " 


30.3 


1.72 


1.62 


90 


Pd(2e) 
Ti(2e) 




0.5 


0.689 
0.201 


0.25 
0.25 


-0.095 




special _B19 (hep) 




1.732 


1.633 


90 


Pd(2e) 
Ti(2e) 




0.5 


0.6666 
0.1666 


0.25 
0.25 






special B19 {B2) 




1.414 


1.414 


90 


Pd(2e) 
Ti(2e) 




0.5 


0.5 
0.0 


0.25 
0.25 




PtTi 


B19 


29.84 


1.758 


1.668 


90 


Pt(2e) 




u.o 


0.6874 
u. lyoo 


0.25 
u.zo 


-0.1512 




B19' 


29.88 


1.762 


1.656 


93.56 


Pt(2e) 
Ti(2e) 


0.0154 
0.4455 


0.6841 
0.1923 


0.25 
0.25 


-0.1561 




Llo 


29.52 


1.000 


1.386 


90 










-0.1397 




bco 


30.57 


1.835 


1.561 


105.80 


Pt(2e) 
Ti(2e) 


0.0907 
0.3502 


0.6816 
0.2001 


0.25 
0.25 


-0.1231 




B19 " 


30.66 


1.75 


1.663 


90 














B19 


29.75 


1.75 


1.666 


90 














B19 " 


30.9 


1.72 


1.62 


90 


Pt(2e) 
Ti(2e) 




0.5 


0.688 
0.197 


0.25 
0.25 


-0.155 



"X-ray, from Ref. 
''X-ray, from Ref. 13. 
'^previous calculation, from Ref. 
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yielding a total of two free internal parameters Wp^/Pt 
and WTi • The unstable mode corresponds to a distortion 
with fixed vp^/pt/vTi ratio. In PdTi, the ratios for the 
phonon and unstable force constant eigenmode as com- 
puted by PWSCF (VASP) are 1.55 (1.48) and 1.75 (1.64), 
respectively. In PtTi, the corresponding ratios are 1.25 
(1.19) and 1.81 (1.62), respectively. 

In Table |, we provide the calculated equilibrium lat- 
tice parameters of the J519 structure, obtained by relax- 
ing all five free structural parameters. The results are in 
good agreement with experiment, aside from the volume 
underestimate typical of LDA. The computed relaxed val- 
ues of the ratio vp^jpijvxi are 1.29 for PdTi and 1.16 
for PtTi. The latter is quite close to the VASP phonon 
result, while the former is significantly lower. To under- 
stand this result more fully, we separate the coupling to 
strain from that to the second (stable) mode of the same 
symmetry by relaxing the structure to _B19 with the lat- 
tice held fixed. The resulting ratios vp^/pt/vTi are 1.51 
for PdTi and 1.41 for PtTi. Thus both couplings are sig- 
nificant in both systems, though the effects fortutiously 
nearly cancel in PtTi. This behavior is compatible with 
continuum models ISltl. 



Comparing the parameter values in Table ^, we see that 
the relaxed B19 structure is close to the ideal hep struc- 
ture that would be obtained from packing monodisperse 
hard spheres, with parameters given in the last line of the 
Table |. The lower symmetry arises from the "decora- 
tion" of the close packed plane with two different atomic 
species, which cannot preserve symmetry and maximise 
unlike near neighbors [Q. This confirms that the phase 
transition is best thought of as a binary equivalent of 
bcc-hcp, not simply as a distortion of the B2 structure. 
There is an interesting contrast with NiTi here: the ex- 
perimentally reported ground state NiTi 519' phase has 
h/a = 1.603 and c/a = 1.417[^, far from hep. 

In the soft-mode approach, there is no guarantee 
that the energy-minimizing freezing-in of one unstable 
mode will stablize the other unstable modes of the high- 
symmetry structure. In the present case, the undistorted 
doubled unit cell contains two sets of (110)(,cc planes each 
of which is unstable to strain-coupled shuffling at all q- 
points. The simplest such mode in the i?19 structure 
is which lowers the symmetry to monoclinic P2i/m. 
The resulting _B19' structure has three additional free 
parameters: the monoclinic angle 7 and two internal pa- 
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FIG. 3: (a) Total energies as a function of uniaxial strain with the internal coordinates being fully relaxed. The lattice 
parameters for A£/esi9=0, 1 are taken from B2 and B19 structures, while others were fixed by interpolation. The "_B2" 
structure here corresponds to freezing in the M-point phonon, and lies 10.5 meV/atom (PdTi) or 30.2 meV/atom (PtTi) 
below the fully symmetric (« = 0) structures. (b)Total energies as a function of monoclinic angle. For B19, B19' and bco all 
degrees of freedom are relaxed. For intermediate angles the lattice parameters a, b, c were fixed by interpolation, while the 
internal coordinates u, v were fully relaxed, u and v vary almost linearly with 7. Symbols indicate calculated values, curve is 
a polynomial fit with dE/d'y = enforced where required by symmetry. 



TABLE II; Optical phonon frequencies at F for PdTi and 
PtTi in the B19 and bco structures by DFPT. The lowest 
frequency mode couples to the strain in the B19-B19' transi- 
tion. 



Alloy 


structure 


Lo (cm ^) 


PdTi 


B19 


73, 90, 101, 139, 159, 184, 199, 231, 251 




bco 


94, 138, 146, 148, 164, 198, 200, 230, 256 


PtTi 


B19 


68, 87, 98, 123, 165, 185, 196, 238, 260 




bco 


119, 171, 172, 173, 185, 210, 225, 247, 284 



TABLE III: Calculated optical phonon frequencies and ra- 
tio of atomic displacements of r4 mode of PdTi and PtTi 
at relaxed _B19 structure from force constant matrix eigen- 
vectors (i.e. ignoring masses). The values in italic are the 
modes leading to 519' structure and are compared with val- 
ues in parentheses which are taken from relaxed 519' struc- 
ture. Frequencies here differ slightly from Table ^ because 
they are calculated from finite displacements using VASP 
rather than DFPT. 



Alloy 


Frequencies 


Upd/Pt/uTi 


PdTi 


74, 103, 159 


-0.280 (-0.217), 3.5684, -1.002 


PtTi 


69, 106, 161 


-0.321 (-0.283), 3.1159, -1.05 



rameters uti and upd- These values are given for the 
relaxed BIQ' structure in Table |[ and compared with 
the normalized eigenvector displacements in Table [IT 
However, it is important to note that the computed T 



phonon modes in the i?19 structures are in fact all sta- 
ble (Table H). The lowering of energy by distortion to 
_B19' cannot be obtained by a pure r4 phonon distortion, 
but only if the strain is allowed to relax simulataneously 
(Table |). This may be the reason that in a previous 
calculation jl^ , _B19 was reported to be the minimum 
energy structure. 

The relative energies of the various relaxed structures 
are given in Table ||. The B\Q total energy is lower than 
B2. A simple estimate of the transition temperature is 
given by /^E = kT^ which suggest of 1050K (PdTi) 
and 1755K (PtTi). These rough values are significantly 
larger than the experimental data for the hysteretical 
transition region[^ ^ ^ (approximately 800K and 
1400K respectively), but show the correct material trend. 
For neither system has a BlW phase yet been observed 
experimentally. The small computed energy differences 
between i319' and B19', translated into temperature, are 
28K and 39K for PdTi and PtTi, respectively This sug- 
gests that the transition to the 519' phase should occur 
at temperatures well below those at which the experi- 
ments were performed, so that our results are fully con- 
sistent with the available experimental work. 

The binary-hep phase interpretation also suggests that 
we should examine the binary- fee equivalent, which is 
the LIq phase. LIq is accessible from B2 by a simple 
(001) shear and has lower energy |Q, however we find 
that LIq has slightly higher energy than 519. We con- 
sider one further structure: at the special 519' values 
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cos 7 — a/26 and 4m — 1 = 2v, orthorhombic symme- 
try is restored. The side of the conventional cell of this 
body centered orthorhombic (bco) (space group Pmna) 
structure is doubled in the b direction, though the prim- 
itive cell still contains four atoms. Although at relatively 
high energy ( |), this provides us a reference point for 
structures with large 7|3^. Note that a further shear to 
cos 7 — a/b would give the _B19 structure once more. 

The potential energies of continuous paths between the 
structures identified above are important for understand- 
ing the transformation mechanism. We compute the en- 
ergies of three paths: B2-B19, B19-S19' and B19'-bco. 
In keeping with the timescale separation between bulk 
strain and atomic motion, we define intermediate con- 
figurations by relaxing the atoms to their minimum en- 
ergy configuration consistent with the applied symmetry 
and strain on the cell. The remaining four strain de- 
grees of freedom are reduced to a single parameter by 
taking interpolations between the strains of the endpoint 
structures. For the B2 structure, we minimize the energy 
assuming the -B19 space group, which gives a smooth evo- 



lution of the structure along the path. /,From Figure 
it is clear that there is no total energy barrier along the 
S2-B19-B19' path, and that B19 represents the total 
energy barrier between i?19' martensitic variants. The 
519 phase can be viewed as the binary equivalent of hep, 
and the transformation path as the binary equivalent of 
the Nishiyama-Wassermann (NW) path. Using the anal- 
ogy with the NW mechanism for the bcc-hcp transition, 
we can attribute the transition to a shuffling of (llO)bcc 
planes 

Structural instabilities in metals are typically related 
to details of the Fermi surface, and we have calculated 
the band structures of i32, B\% and B\%' to investigate 
this. In PdTi/PtTi the band structure is dominated by 
the rf-bands, with the Pd/Pt bands lying below the Fermi 
level and almost fully occupied, and the Ti d-bands lying 
above the Fermi level (Figure ^ and ^), the band cen- 
ters being offset by some 6eV. The free-electron like sp 
bands are very broad, and play little role in the bonding 
except to donate some electrons to the Ti-d band. The 
large strain involved in the B2 — _B19 transition means 




that the Fermi surfaces are quite different {B19 having 
the lowest DOS at iJi?]!^) and this transition cannot be 
related directly to the band structure. By contrast, the 
B19-B19' transition is accompanied by the opening of a 
pseudogap at the Fermi level, a typical signature of in- 
creased stability. The band structures are very similar, 
the small difference which stabilises B19' being traceable 
to the shifting above the Fermi level of a pocket of elec- 
trons around B. The distortion to B19' is just enough 
to complete the topological phase transition which elim- 
inates this pocket of electrons in both materials. 



IV. CONCLUSIONS 

In conclusion, we have performed ab initio calculations 
of the structural energetics of PdTi and PtTi. In each 
case we predict that the low temperature ground state 
structure will be B19', with the (observed) B2 and B19 
phases being dynamically stabilized. There are no total 
energy barriers between the structures, meaning that the 



phase space microstates that belong to the B19' structure 
also belong to the B19 and B2 structures. 

In contrast to NiTi|jl^ , the entire T—M phonon branch 
is unstable. We showed that the B19 structure can be ob- 
tained by a "freezing in" of phonons of the B2 structure 
coupled to the shear associated with the c' = (cn — ci2)/2 
elastic constant, but that no single dynamical-matrix 
or force-constant-matrix eigenvector leads to the low- 
symmetry phase. 

The B19' then corresponds to a further strain coupled 
to a r4 phonon of the B19 phase. Tracing the atomic 
motions of these instabilities shows that they are both 
related to shears of alternate (110)b2 phases, and hence 
that the transition mechanism is the binary equivalent of 
the Nishiyama-Wassermann bcc-hcp mechanism. 
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